New facts about muon production in Extended Air Shower simulations 
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Whereas air shower simulations are very valuable tools for interpreting cosmic ray data, there is a 
long standing problem: is seems to be impossible to accommodate at the same time the longitudinal 
development of air showers and the number of muons measured at ground. Using a new hadronic 
interaction model (EPOS) in air shower simulations produces considerably more muons, in agreement 
with results from the HiRes-MIA experiment. We find that this is mainly due to a better description 
of baryon-antibaryon production in hadronic interactions. This is a new aspect of air shower physics 
which has never been considered so far. 
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Since more than ten years detailed extended air shower 
(EAS) simulations play a decisive role in interpreting 
measurements from ground based cosmic ray measure- 
ments. This concerns for example the chemical composi- 
tion of cosmic rays in the KASCADE experiment [ll or 
the primary energy determination in AG AS A array 3|- 

An air shower is initiated by a very energetic proton 
or nucleus (primary particle), which interacts with air 
by producing many secondary hadrons, which interact 
again, and so on. Neutral pions play a special role, since 
they decay into gammas which initiate an electromag- 
netic shower each. The latter one is well under control 
(elementary processes of QED), whereas the hadronic in- 
teractions require models, being tested against accelera- 
tor data - at much lower energies though than the high- 
est primary energies. It turns out that more than 90% 
of the initial energy goes into the (well-known) electro- 
magnetic part, whereas the rest shows up as muons from 
the hadronic decays. The muonic part depends therefore 
strongly on the hadronic modeling (up to a factor of 2 
difference in different models) , while the longitudinal de- 
velopment of the electromagnetic part (and in particular 
its maximum Xmax) are relatively robust (less than 10% 
variations between models). 

Using the currently employed hadronic interaction 
models (QGSJETOl [|, QGSJET II [1], SIBYLL 2.1 0]), 
one has a tendency to have less muons in the simulations 
than observed by the experiments. Direct measurements 
of high energy muons (100 TeV) [1, 0, Bll as well as 
experiments using low energy muons (1 GeV) like KAS- 
CADE [H or HiRes-MIA [lO|, show inconsistencies be- 
tween experimental data and simulations. Furthermore, 
at very high energy, the Pierre Auger Observatory finds 
a discrepancy between the energy reconstruction of the 
primary cosmic rays using a purely experimental method 
or using a method based partly on the muon density at 
ground [ll| and air shower simulations. The 25% differ- 
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ence could be explained by a lack of muons in the sim- 
ulations. Many attempts have been made to force the 
models to increase the muon production without chang- 
ing Xmax (well constrained by data) without success [ij] ■ 

It should be noted that the number of muons predicted 
from shower simulations has important consequences on 
the astrophysical interpretation of the very high energy 
cosmic ray spectrum. If the muon number currently 
predicted is right, "new physic" has to be introduced 
to explain the observation of events beyond the "GZK- 
cutoff" dl of lO^OeV, observed by the AGASA experi- 
ment (for a review see [T^|). 

In this work, we discuss the consequences of introduc- 
ing the recently developed high energy hadronic interac- 
tion model EPOS into the air shower simulation models 
CORSIKA :1&] and CONEX [ll,[l3,[il. EPOS is a con- 
sistent quantum mechanical multiple scattering approach 
based on partons and strings [l9|, where cross sections 
and the particle production are calculated consistently, 
taking into account energy conservation in both cases 
(unlike other models where energy conservation is not 
considered for cross section calculations [23|). A special 
feature is the explicit treatment of projectile and target 
remnants, leading to a very good description of baryon 
and antibaryon production as measured in proton-proton 
collisions at 158 GeV at CERN [21]. Motivated by the 
very nice data obtained by the RHIC experiments, nu- 
clear effects related to CRONIN transverse momentum 
broadening |22j . parton saturation, and screening have 
been introduced into EPOS [1^. Furthermore high den- 
sity effects leading to collective behavior in heavy ion col- 
lisions are also taken into account It appears that 
EPOS does very weU compared to RHIC data [H, Hi], 
and also all other available data from high energy particle 
physic experiments (ISR,CDF and especially SPS exper- 
iments at CERN) [27|. As a result, it is the only model 
used both for EAS simulations and accelerator physics 
which is able to reproduce consistently almost all data 
from 100 GeV lab to 1.8 TeV center of mass energy, in- 
cluding antibaryons, multi-strange particles, ratios and 
pt distributions. It is intensively tested against data, but 
as explained, it is not a simple fit of data. And in particu- 
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lar, since this model is used for accelerator physics, many 
data are used to constrain the model parameters which 
are not a priori linked to cosmic rays and air showers. 

For our analysis, CONEX and CORSIKA are 
used to simulate the air shower development, using 
GHEISHA [2^ as low energy hadronic interaction model 
below 80 GeV. For the high energy interactions (above 
80 GeV) EPOS 1.35 is used, and as a reference the most 
commonly used interaction model QGSJETOl. 

One of the most important observables in air shower 
physics is the distribution of the electron number as 
a function of the depth X, the latter one representing 
the amount of air traversed by the shower, expressed in 
g/cm^. The maximum X^ax of this distribution is a func- 
tion of the energy of the primary particle, as shown in 
Fig. [H where the results represent averages over many 
showers. The experimental data from HiRes [l^ (points 
in Fig. [T|) refer to unknown primary particles, therefore 
one usually compares the data with the two extremes 
(protons and iron) from simulations. Here we show re- 
sults for EPOS (full lines) and QGSJETOl (dotted lines). 
The upper lines represent protons and the lower lines 
iron. Both models are compatible with experimental 
data which seems to show a lightening of the primary 
cosmic ray composition between 10^^ eV and 10^^ eV. 
But even at the lower energies, the average primary par- 
ticle doesn't seem to be heavier than a carbon nucleus 
as shown by the dashed line calculated with EPOS for 
carbon-induced showers. Above 10"'^^ eV, data are well 
reproduced by proton induced shower. 

A complementary observable is the muon number at 
ground, for example expressed via the density /9^(600) 
of muons per squared meter at a lateral distance of 600 
m from the shower core (impact point) as measured by 
the MIA detector as a function of the primary en- 
ergy, shown in Fig. [2] as triangles. Again we show 
shower simulations for protons and iron for EPOS and 
QGSJETOl, using the same conventions as in Fig.[TJ The 
HiRes-MIA data are now compatible with the EPOS re- 
sults, using a heavier primary composition (but not too 
heavy, more like carbon) at 10^^ eV and a lighter one 
(proton) at 10^^ eV. Compared to QGSJETOl, it is a 
shift of about 40% in the number of muons at ground. 
The proton line from EPOS lies exactly on the iron line 
from QGSJETOl (which was the model giving the highest 
number of muons before). 

So for the first time, both Xmax and muon data are 
compatible with a change of the average incident particle 
from carbon to proton, between 10"'^'^ eV and 10^® eV. 

Not only the absolute value of the muon density has 
changed but also the slope is slightly higher (less nega- 
tive) in EPOS compared to QGSJETOl. This can be seen 
on a larger energy scale Fig. O the number of muons ar- 
riving at ground divided by the primary energy in GeV is 
shown as a function of the primary energy between 10^"* 
eV and 10^^ eV. The EPOS curves are much flatter: at 
the lowest energy, the EPOS proton line is at most 25% 
higher than QGSJET, but at 10^° eV EPOS is a fac- 
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FIG. 1: The shower maximum Xmax as measured by the 
HiRes-MIA 10] and the HiRes collaboration 29] as a function 
of primary energy, compared to proton and iron induced show- 
ers simulated with EPOS (full lines) and QGSJETOl (dotted 
lines) as high energy hadronic interaction model. The thin 
line is calculated with EPOS for carbon nuclei as primaries. 
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FIG. 2: The muon density p^(600) from the MIA exper- 
iment [l(|] (triangles) as a function of the primary energy, 
compared to simulated proton and iron induced showers with 
EPOS (full line) and QGSJETOl (dotted line) as high en- 
ergy hadronic interaction model. For better comparison with 
Fig. 1, the dashed line correspond to carbon induced showers 
with EPOS. 



tor 2 higher and gives even more muons with a primary 
proton than QGSJETOl for iron induced showers. As a 
result, EPOS should be compatible with the KASGADE 
data. It will most probably give a lighter average compo- 
sition, but a precise answer has to await full simulations 
(in preparation). However, EPOS will give very different 
results for the AUGER spectrum based on surface de- 
tector for instance compared to QGSJETOl simulations. 
In addition, EPOS has a slightly higher elongation rate 
(slope of X,nax as a function of primary energy) compared 
to QGSJETOl but stays fully compatible with the HiRes 
data at highest energies. 
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FIG. 3: Total number of muons at ground divided by the 
primary energy expressed in GeV as a function of the primary 
energy for proton and iron induced shower using QGSJETOl 
(dotted lines) or EPOS (full lines) as high energy interaction 
model. 



For a deeper understanding of what makes the differ- 
ence between EPOS and older models, one has to have 
a closer look at the EAS physics. Looking at all parti- 
cles produced in EAS by EPOS, we find that not only 
the number of muons is increased but also the number of 
baryons and antibaryons - referred to as (anti)baryons 
in the following. As a test, we modified EPOS pa- 
rameters artificially to reduce the number of produced 
(anti)baryons by a factor of two, hence not reproducing 
accelerator data anymore. As a result, the number of 
muons is reduced on the average by 30%, being quite 
close to QGSJETOl. At the same time, Xmax is not 
changed. So we find - and this has never been discussed 
in the literature so far - that (anti)baryon production is 
a very efficient mechanism to affect the muon numbers 
without touching Xmax- 

As a cross check, we increased artificially the 
(anti)baryon production in the SIBYLL 30] hadronic in- 
teraction mode (commonly used for EAS simulations) by 
a factor of 2, resulting in 30% more muons. 

Up to now, the correlation between the number of 
(anti)baryons (about 1% of the particles in the hadronic 
shower) and the number of muons in air showers has 
never been studied. But, thanks to a simple Heitler 
model generalized to hadronic showers jsT], [34I , one can 
easily understand the role of the antibaryons in EAS. 
In this kind of toy model, a hadronic interaction of a 
charged particle with energy E will produce iVtot new 
particles with energy E/Ntot, with TVem particles (mainly 
TT^) transferring their energy to the electromagnetic chan- 
nel via immediate decay into photons, each one initiating 
an EM shower. The other A^had = A'tot — A'em parti- 
cles re-interact with air after a constant interaction path 
length Ahad and thus contribute to the hadronic cascade. 
In the electromagnetic shower, exactly 2 secondaries with 
equally shared energy are produced at each interaction. 



after a constant interaction path length Aem- Introduc- 
ing a characteristic energy {E^ec = 150 GeV), where pi- 
ous are assumed to decay into muons, the number of 
muons for a shower with primary energy Eq can be writ- 
ten as [33|] 



= {A^had}" 



Eq 
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with a = In A'had/ In A'tot < 1, and where n is the number 
of hadronic generations in the shower. Introducing R = 

A^had/Atot, we have 
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Eq. ^ shows that the muon number depends strongly on 
the ratio R of the number hadrons initiating a hadronic 
sub-cascade to the total multiplicity, which is under- 
standable since the difference between these two quan- 
tities are particles giving all their energy to the electro- 
magnetic channel - not producing muons. 

Usually these kind of toy models consider only pions as 
secondary particles. As a consequence, one has R — 2/3. 
So the muon number depends only on A'tot- The latter 
one affects as well Xmax, since from our toy model we 
obtain [331 



Ahad + A 



EM 



In 



^0 



Atot-Bc 
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where Ec — 85 MeV is the critical energy (energy of the 
electromagnetic particles at the shower maximum in air 
and energy of particles disappearing from the shower). 
So EAS simulations based on two different interaction 
models, producing different total multiplicities, should 
disagree for both Amax and muon numbers. 

Let us now be more realistic, and consider all kinds 
of hadrons, including (anti)baryons. Particle production 
in hadronic interactions is model dependent, and so is 
the precise value of R. With R being less than 1 and 
A'tot >> 1, the muon number depends very sensitively 
on the ratio R. One may imagine two interaction mod- 
els with the same Atot but with more baryons in one 
model compared to the other, corresponding to smaller 
A^EM and bigger R values. With R being even slightly 
higher, a is closer to one, increasing both the number 
of muons and the slope as a function of the energy as 
observed Fig. [31 This muon increase is also intuitively 
understandable: more baryons in one model compared 
to the other lead to a smaller Aem, less energy is trans- 
fcrcd to the electromagnetic component. There are more 
hadronic generations and thus more muons are produced. 

As a result, it is clear that variations of the number of 
baryons and antibaryons affect i?, but why is the effect 
so strong although baryons and antibaryons represent 
only 1% of all hadrons. One should not forget that only 
hadrons with a large longitudinal momentum fraction x e 
(>0.1) contribute significantly to the cascade. Particle 
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FIG. 4: EPOS 7r° energy spectra for ir-Air interactions (full 
line) and for p-Air interactions (dashed line) at IC GeV ki- 
netic energy. 

yields at large x e are quite different from those at central 
rapidities, measured in accelerator experiments. Fig. [J] 
shows the energy spectra of 7r° produced in p-Air (dashed 
line) and n-Air (full line) reactions at 10^ GeV kinetic 
energy for EPOS, tt'^'s produced in p-Air interactions are 
much softer than in T:-Air simply because p-Air interac- 
tions do not produce leading tt" whereas n-Air do so 
because of charge exchange. As a consequence, the ratio 
R « 1 — A^TfO /Ntot for particles with more than 50% of the 
interaction energy is close to 15% for ir-Air but about 1% 
for p-Air. In other words, producing only slightly more 
baryons per hadronic interaction in one model compared 
to another will decrease considerably the number of fast 
neutral pions (w iVEivi) in the next generation, leaving 
more contributors to the hadronic cascade and providing 
thus more muons. It is worthwhile to say that having by 
any means more leading 7r°, will decrease the number of 
muons and the slope as a function of the energy [34] . 
Why does the muon yield in EPOS differ so much from 



other models? As already said, EPOS has not been de- 
signed to be used for EAS simulations, but to understand 
hadronic interactions. Therefore all available data on p- 
p, p-A, A- A, TT-p/A, K-p/A have been considered, in par- 
ticular data on observables like (anti)baryons which so far 
have not been considered to be important for EAS simu- 
lations. Differences between hadronic interaction models 
used for EAS simulation will be discussed in detail in a 
forthcoming publication. 

To summarize: simulating air showers by using the 
new high energy hadronic interaction model EPOS re- 
sults in an increase of the muon density at 10^^ eV of 
about 40% compared to QGSJETOl calculations. So for 
the first time, both Xmax and muon data are compatible 
with a change of the average incident particle from heavy 
to light element, in the energy range between 10"'^^ eV and 
10^* eV. It was shown for the first time that (anti)baryon 
production plays a much more important role in EAS 
physics than expected. An increased (anti)baryon pro- 
duction (in one model compared to another) increases the 
number of interactions where no leading 7r° is produced, 
more energy goes into hadronic sub-showers, leading to 
more hadron generations, and finally to more muons. 

The presented EPOS results have certainly to be con- 
firmed by comparing to other air shower experiments 
(work in progress), but it seems already now that the 
Cosmic Ray energy spectrum based on EAS simulation 
will have to rescale its energy to lower values. Another 
consequence: since there is more energy in the hadronic 
part of the shower, the calorimetric energy measured by 
fluorescence detectors such as HiRes will correspond to a 
larger primary energy. The conversion factor is changed 
by about 5%. A more quantitative statement on the CR 
spectrum requires intensive simulations which will be pre- 
sented soon. 
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